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' Abstract. Hartree-Fock-Bogolyubov (HFB) calculations making use of a recently proposed microscopic effective pairing 

interaction are presented. The interaction was shown to reproduce the pairing properties provided by the realistic AVI 8 force 
very accurately in infinite matter. Although finite-ranged and non-local, it makes 3D HFB calculations in coordinate space 
tractable. As a first application, basic pairing properties of calcium isotopes in their ground-state are studied. By comparing 
the results with those obtained using a standard Density-Dependent Delta Interaction, the crucial isovector character of the 
microscopic interaction is highlighted. 

in ' 

1. INTRODUCTION 

CO ' 

The structure of the nucleus and the properties of extended nuclear systems strongly depend on their possible superfiuid 
nature. In finite nuclei, pairing constitutes the main part of the residual interaction and has a strong influence on most 
of low-energy properties of the system [1]. In extended systems such as neutron stars, pairing is a decisive ingredient 
of dynamical and thermal evolutions [2, 3]. Despite its major role, the present knowledge of the pairing force and the 
nature of pairing correlations in nuclei, that is, the way Cooper pairs are formed in the nuclear medium out of the 
""^ ' strong nucleon-nucleon (AW) force, is quite poor. Properties such as the range of the effective pairing interaction, its 

link to the bare force, its possible surface character in finite nuclei and its density dependence (in particular isovector) 
still have to be clarified [4, 5, 6, 7]. 

Regarding self-consistent mean-field calculations of finite nuclei, only phenomenological forces such as the Gogny 
force [8] or (Density-Dependent) Delta Interactions ((DD)DI) [5, 6, 9, 10] have been used in the pairing channel so 
far. One exception exists however [11], where the Bardeen-Cooper-Schrieffer (BCS) gap equation was solved in the 
Hartree-Fock (HF) basis of 120 Sn using realistic AV 18 bare force [12]. It required the treatment of single-particle states 
up to 800 MeV! Although successful in describing low-energy nuclear structure over the known mass table [1], the 
Gogny force and DDDI lack a clear link to the bare nucleon-nucleon (AW) interaction 2 . This feature strongly limits the 
reliability of their analytical structure such as their possible density dependence. Also, their direct fit to nuclear data 
through mean-field calculations makes probable the re-normalization of beyond-mean-field effects. This is a significant 
limitation if one wants to go explicitly beyond that level of approximation. Finally, their fits performed onto very 
limited sets of nuclei around stability make their extrapolated use toward the drip-lines unsafe. For instance, while 
such phenomenological forces all provide similar and reasonable pairing properties around stability, the predicted 
location of the two-neutron drip-line can differ by up to ten/twenty mass units depending on the force used [16]. 

To improve on that situation, a microscopic effective interaction explicitly linked to the bare AW force, and 
equivalent to it at the mean-field level, was proposed recently to treat pairing correlations in the l So channel [15]. 
BCS pairing properties provided in infinite matter by AV 18 were reproduced very accurately. These properties dealt 
not only with the gap at the Fermi surface as a function of density, but also with the momentum dependence of the gap 
at fixed density [15]. 

In the present paper, we discuss the first results of 3D HFB calculations using that force and we focus on the isovector 
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- Even if some significant differences remains, one can argue that the Gogny force behaves almost like a bare force in the So and ' Di channels, 
especially when the D1S [13] parameterization is used [14, 15]. 



and low-density properties of the interaction. In section 2, we briefly outline the formalism and the characteristics of 
the new pairing force. Results obtained along the calcium isotopic chain are discussed in section 3. They are compared 
to those obtained using a standard surface peaked DDDI and a zero-range approximation to the new microscopic force. 
Conclusions are given in section 4. 



2. FORMALISM 

To treat pairing, one needs to specify the many-body technique used and the appropriate interaction to insert into 
the calculation at the chosen level of approximation. The latter depends both on the situation and on the system. 
In the present case, we concentrate on a (self-consistent) mean-field description of finite nuclei using the HFB 
method. Eventually, calculations beyond the mean-field are to be performed. Typically, correlations associated with 
symmetry restorations (Projected Mean Field Method) and large amplitude motion (Generator Coordinate Method) 
are considered [17, 18, 19]. Thus, one has to identify the appropriate vertices to be used coherently at each level of 
approximation. 

While variational calculations are of no direct help in that respect, perturbative methods using Green's function 
or Goldstone formalisms provide guides to do so. In particular, such many-body theories show unambiguously that 
the interaction to be used in the particle-particle (p-p) channel at lowest order in irreducible vertices is the bare AW 
force 3 . At the next order, the irreducible pairing vertex involves the so-called polarization diagrams. This situation is in 
contrast to the particle-hole (p-h) channel where one may need to regularize the repulsive core of the bare interaction 
from the outset through the definition of an in-medium two-body vertex like the G-matrix [20]. This stresses the fact 
that the effective forces may differ in the two channels at a given level of approximation. 

The approximate ground-state energy E HFB of a nucleus is a functional of the one-body density matrix p?. = 

( <£> | cj Cj | <t>) and pairing tensor K% = ( <£> | c/ Cj | <E>), where | <I>) is the HFB state 4 . For a general presentation of the HFB 
formalism, we refer to Ref. [17]. Also, the two-basis method used to solve the HFB problem iteratively is discussed in 
Ref. [21]. Finally, a detailed presentation of the method applied to the new interaction will be soon available [22]. 



The pairing field A?, reads as: 



aK ij z kl 

where, as already explained, V >s ° is the bare NN force. In Ref. [15], the effective vertex was then introduced by 
recasting the previous gap equation written in the canonical basis 5 into a fully equivalent expression [15]: 
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where ^'^(O) is an off-shell in-medium two-body matrix summing p-p and h-h ladders in the superfluid phase. 
In Eq. 2, the expression of the pairing tensor in its canonical form jc mm = A^/2E m was used, where E„ = = 



[hm — Hi) 2 + Am 2 , while h q m is the diagonal matrix element of the HF field in the canonical basis and ji q the chemical 
potential. Identifying the two sides of Eq. 2, the effective pairing vertex is naturally defined through its antisymmetrized 
matrix elements in the canonical basis of the Bogolyubov transformation as: 



(V$) =2p«(^>)),.- (3) 
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3 This constitutes our motivation for the mean-field level. With the G-matrix [20] in the p-h channel, such a mean-field theory precisely aims at 
treating the nucleus as a system of independent pairs, including the correlations associated with the existence of Cooper bound-states in the medium. 

4 The quantities are defined in an arbitrary basis and the isospin quantum number q (1/2 for neutrons and -1/2 for protons) is specified to make 
clear that isospin mixing is considered neither in the p-h channel nor in the p-p channel. 

5 The canonical basis corresponds to the single-particle basis diagonalizing the one-body density matrix p? = p? <5 ; , = pj and putting the pairing 
tensor in its canonical form (c? ; = (A <5,j. pO, plays in the canonical basis the role of the usual BCS occupation number. Which states are paired 
is a by-product of the Bogolyubov transformation solution of the problem. States are not paired a priori as in the BCS approximation. 
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In Ref. [15], ^ ° (0) was calculated explicitly in infinite matter starting from a separable form of the bare AW 
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interaction in the 1 Sq channel. It was studied in detail and shown to take a closed form in coordinate space: 

,A \-P a f e-Eii \7-7i\ 2 /2a 2 

9 s o (0) (n, r 2 , r 3 , r 4 ) = X s ° J dr f[p q {r)\ 

where P a is the spin-exchange operator 6 , while a = V0.52 fm and A' s ° = - 840 MeV.fm 3 denote the range and 
the intensity of the force, respectively. No further adjustment is to be made in finite nuclei. The functional f[Pq{f)\ 
incorporates the density dependence of the T z = 2q component of the effective interaction [15]. The density dependence 
stems from the re-summation of p-p and h-h ladders in the medium. The effective vertex is thus finite ranged, non 
local, total-momentum dependent and density dependent 7 . However, the computing cost of corresponding 3D HFB 
calculations is, through the two-basis method [21, 23], of the same order as for a zero-range interaction. 

By re-summing the effect of pairs scattered at high-energy into the effective vertex, the latter is soft even if the bare 
force has a hard core. Also, a smooth cut-off 2 pf n emerges naturally in the gap equation through its recast. This cut-off 
further limits the necessity to use large basis sets as in Ref. [11] and makes zero-range approximations of the effective 
vertex meaningful. The pairing problem is regularized in a similar way to what was proposed in Refs. [7, 24, 25] 
using re-normalization techniques. It is worth noting that the derived cut-off differs from all ad-hoc ones used in 
connection with usual DDDI. Thus, a zero-range (ZR) approximation providing identical gaps at the Fermi surface in 
infinite matter was defined in Ref. [15]. The coefficients entering the functional f[pq(j)\ differ from the ones used in 
the finite range case. Performing such an approximation, the roles of the range and of the density dependence of the 
interaction could be disentangled [15]. In particular, the surface-enhanced character of phenomenologically optimized 
DDDI [5, 6] was demonstrated and shown to be, to a large extent, a way of re-normalizing the range of the interaction. 
Also it was shown that usual DDDI miss the low-density behavior of the effective pairing force. 

Both the finite-range microscopic force and its zero-range approximation depend on the density p q (j) of the 
interacting particles rather than on the total matter density. It is to be contrasted with usual DDDI which often take the 
form: 
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<S(n - r 3 ) 8(r 2 - r 4 ) <S(n - n) , (5) 



Pc 

where po(r) denotes the total matter density (the local scalar-isoscalar part of the one -body density-matrix) while p c 
is equal to (one-half) the saturation density for the surface (half-surface) type pairing force. 

The role of the finite range, the isovector density-dependence and the low density behavior of the pairing force, as 
well as the regularization scheme used together with contact approximations has to be addressed in detail. While we 
briefly focus on the isovector properties of the interaction in the present communication, we refer to Ref. [22] for an 
extensive study of all other issues. 



3. RESULTS 

We performed 3D HFB calculations of Ca ground-states from proton to neutron drip-lines [21]. Good particle-numbers 
were approximately restored before variation through the Lipkin-Nogami (LN) procedure [26]. Self-consistent block- 
ing and time reversal symmetry breaking were included in the calculations of odd isotopes. The Sly4 Skyrme force [27] 
was used in the p-h channel. Each calculation was repeated three times using the microscopic finite range force defined 
through Eq. 4, its zero range approximation and a standard surface peaked DDDI as given by Eq. 5 [10]. The latter 
was adjusted together with a phenomenological cut-off defining an active window (±5MeV) around the Fermi energy. 

By keeping the force in the p-h channel fixed, we probe the vertex used in the pairing channel, including the self- 
consistent coupling between the two channels. Of course, properties of the force in the p-h channel have an impact 
on the results. In that respect, it is worth noting that the considered DDDI was adjusted on properties of (non-exotic) 
nuclei using the SLy4 parameterization in the p-h channel [10], and thus, is consistent with the isoscalar effective mass 



The force acting only in the relative S-wave, the projection on the spin singlet corresponds to a simultaneous projection on the isospin triplet. The 
T = 1 neutron-proton pairing (T z = component) is not considered here. 

7 Although the effective vertex breaks total-momentum and angular-momentum conservation of the interacting pair, as well as isospin symmetry, it 
does so in such a way that the energy functional itself remains invariant under rotation in isospin and real spaces, as well as under translation [22]. 




FIGURE 1. Left panel: one-neutron separation energy Sn(A) = E(N — 1,Z) — E(N,Z) in Ca isotopes. Right Panel: two-neutron 
separation energy S2n(A) = E(N — 2,Z) — E(N ,Z) for the same nuclei. Experimental [29] (diamonds) and theoretical values for 
three different pairing forces are displayed. 



m* jm = 0.7 predicted by the latter. The microscopic pairing forces were adjusted once for all without any reference 
to finite nuclei with the property of not depending explicitly on the effective mass appearing in the p-h channel [15]. 
However, according to our definition of the mean-field, they should be used at the HFB level together with a p-h vertex 
providing an effective mass consistent with a G-matrix supplemented by a three-body force [28]. 

In Fig 1, calculated one and two-neutron separation energies are compared to experiment [29]. Overall, the agree- 
ment is (at least) of the same quality as the one obtained from other calculations of the same type [7, 30]. There are 
interesting differences which are beyond the scope of the present discussion. We rather concentrate on the effect of 
the pairing force used in the region of unknown exotic nuclei, because it is where its choice is of crucial importance. 
As an example, its influence on the position of the two-neutron drip-line is clearly seen on the right panel of Fig. 1 . 
While the usual DDDI predicts it to be located at N — 44, the stability against two-neutron emission extends up to 
N = 50 when the microscopic forces are used. Note that a difference of two mass units in the predicted position of the 
drip-line for such light semi-magic nuclei can translate into a difference of ten mass units for lead isotopes [16]. 

Average neutron and proton pairing gaps calculated with the three pairing forces are plotted in Fig 2 along the Ca 
isotopic chain 8 . We see non-trivial differences between the predictions of the phenomenological DDDI and of the 
microscopic forces. We refer to Ref. [22] for a detailed discussion. Here we simply stress the different isovector trend 
of those predictions. While the magnitude of the neutron gaps are of the same order for nuclei around the stability line, 
the phenomenological DDDI provides much too strong pairing in neutron rich nuclei as anticipated in Ref. [15]. The 
overshoot of simple DI by usual DDDI near the drip-line was also identified [16] and often thought to be a result of 
the surface-peaked character of the latter 9 . In fact, in both cases, the primary cause of the overshoot is the improper 
isovector character of usual surface-enhanced DDDI associated with their dependence on the total (isoscalar) density 10 . 
Indeed, such a DDDI, when adjusted on nuclei having very similar neutron and proton densities, will provide stronger 
(weaker) pairing in a region of neutron rich matter than a force independent of the density or depending on the neutron 
(proton) density. This situation is also illustrated in the right panel of Fig. 2. While the LN prescription is responsible 
for the non-zero value of the proton gap at Z = 20, the latter should not evolve much with neutron number. However, 
< A P (A) > presents an artificial slope when using the DDDI whose origin goes precisely along the line of the previous 
argument. When repeating the calculation without LN, the proton gap even switches on artificially for N < 24 in the 



Values appear for odd and even particle numbers. For odd particle numbers, no blocking was considered here as such HFB states constitute the 
proper reference on top of which the blocking as to be eventually performed to describe the final odd state [31]. The average gaps are not compared 
with experiment because we do not want to discuss here how they are related to odd-even mass differences [32]. 

9 The way each interaction couples to the continuum also plays a role. Usually, pairing correlations are diminished by a strong coupling to the 
continuum. The way the continuum is treated numerically also influences the results near the drip-line [33]. 

10 It is proved here not to be related to the range of the interaction. The gaps obtained with the finite range force do not differ significantly from 
those obtained with its zero-range approximation. 
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FIGURE 2. Neutron (left panel) and proton (right panel ) average pairing gaps < A q (A) >= A„ K% n / K% n calculated along 
the Ca isotopic chain with three different pairing forces. 
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FIGURE 3. Neutron pairing gap at the (neutron) Fermi surface as a function of the neutron density. Left panel: symmetric nuclear 
matter. Right panel: pure neutron matter. The results obtained with the three pairing forces used in the present paper are compared 
to the results derived from the realistic AVI 8 AW force [12]. To make the theoretical comparison clear, all gaps are calculated with 
free kinetic energies as single-particle energies. 



case of the DDDI. 

In order to confirm that interpretation, the BCS neutron gap at the Fermi surface in symmetric nuclear matter and in 
pure neutron matter is plotted in Fig. 3 as a function of the neutron Fermi momentum k n F . Because we use free kinetic 
energies as single -particle energies, A"(kp) obtained from AVI 8 or from the density-dependent microscopic forces are 
the same in symmetric matter and in pure neutron matter. On the other end, while the DDDI adjusted on stable nuclei 
reproduces rather well the gap in symmetric matter, it strongly overshoots it in neutron rich matter. This is due to its 
dependence on the total density and cannot be related here to any surface effect. Furthemore, beyond mean-field effects 
in the bulk one could eventually incorporate will always lead to a decrease of A"(kp) in infinite neutron matter [34]. 

As a last point, we discuss the low-density behavior of the pairing force. It was shown in Ref. [15] that the intensity 
of the effective microscopic force strongly rises at low density due to the very large scattering length of the AW 
interaction in the 1 Sq channel. The influence of this strong attraction can be seen in the right panel of Fig. 1 where 
nuclei are stabilized beyond the sub-shell N — 40 and the two-neutron drip-line pushed back by ten mass units. Indeed, 
the gaps obtained from the microscopic force resist their otherwise decreasing trend near the drip-line because of the 
increasing importance of low densities. Such an effect is not seen for the (still too strong) gaps calculated from the 
DDDI. Strong low-density dependence of DDDI were simulated phenomenologically in Ref. [6] and shown to bring 



about pathologies. This highlights the fact that such a density dependence should be used in connection with the 
corresponding microscopically derived cut-off. 

4. CONCLUSIONS 

We presented results of the first (3D) HFB calculations performed in finite nuclei using a recently proposed micro- 
scopic effective pairing interaction. The isoscalar and isovector density-dependences derived ab-initio provide the 
pairing force with a strong predictive power when extrapolated toward the drip-lines. We concentrated here on that 
aspect by studying basic pairing properties of calcium isotopes in their ground-state. By comparing the results with 
those obtained from a standard Density-Dependent Delta Interaction, the crucial isovector character of the microscopic 
interaction was highlighted. 

In the near future, more systematic HFB calculations will be presented to identify the role of the bare AW force in 
building pairing in finite nuclei. Through comparisons with experiment, an indirect measure of missing correlations in 
the p-p channel will be realized. Local theories of pairing will be challenged by probing the importance of the finite 
range of the force, especially when describing low-energy excited states. Later, the use of the microscopic force will 
be extended to beyond mean-field calculations. 
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